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^ ■ Abstract 



The inclusive cross section for the production of J/ip in e + e annihilation is 



Ph. studied using a recently developed factorization formalism. In addition to the 



conventional color-singlet contribution, we include contributions from cc pairs 



^ . that are produced at short distances in color-octet states. The color-octet 

contributions are suppressed by the smaller probability for the color-octet cc 
pair to evolve into a state that includes a i/j, but they may dominate near the 
upper endpoint of the ip energy spectrum. The signature for the color-octet 
contribution is a dramatic change in the angular distribution of the tp near 
the endpoint. 
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Charmonium is a bound state of a charm quark and antiquark held together by the 
strong force of quantum chromodynamics (QCD). Because QCD confines color charges, the 
quark and antiquark must be in a singlet state with respect to color. The production of 
charmonium requires the creation of a cc pair with energy greater than 2m c , where m c is 
the mass of the charm quark. Since the QCD coupling constant is relatively weak at the 
scale 2m c , the production of charmonium should be amenable to theoretical analysis using 
perturbation theory in a s . Most theoretical predictions of charmonium production have been 
based on the "color- singlet model" [l]. In this model, the cc pair that forms charmonium 
is assumed to be produced in a color-singlet state by parton collisions whose cross sections 
can be calculated using perturbation theory. The probability that a cc pair produced in a 
color-octet state will bind to form charmonium is assumed to be negligible. 

Recent developments in both theory and experiment have challenged this simple assump- 
tion. Bodwin, Braaten and Lepage have developed a factorization formalism that allows 
the systematic calculation of inclusive cross sections for charmonium states to any order in 
the QCD coupling constant a s and to any order in v 2 , where v is the typical relative velocity 
of the charm quark. This formalism implies that "color-octet processes" , in which the cc pair 
that forms charmonium is produced at short distances in a color-octet state, must contribute 
to the cross section at some level. The assumption that charmonium production is always 
dominated by color-singlet processes has also been challenged by recent experimental data. 
The CDF collaboration has measured the cross sections for prompt charmonium production 
at large transverse momentum at the Tevatron || , using a silicon vertex detector to separate 
the prompt signal from the background of charmonium produced by the decay of B hadrons. 
At the large values of pt available at the Tevatron, the dominant production mechanism is 
fragmentation |4[] , the formation of charmonium within the jet of a high energy parton. For 
the J PC = 1 charmonium states J/ip and ip', the prompt cross sections were measured to 
be more than an order of magnitude larger than the predictions of the color-singlet model, 
even after including the effects of fragmentation [Q. 

One of the proposed explanations for the CDF data on prompt ip production is that the 
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large cross section comes from a color-octet term in the fragmentation function for a gluon to 



short distances of order l/m c or smaller, with the binding of the cc pair into a color-singlet 
ip occurring at longer distances. This explanation reproduces the pr-dependence of the CDF 
data. The normalization of the data can be fit by adjusting the value of a matrix element 
that measures the probability for the formation of ip from a pointlike cc pair in a color-octet 
Si state 0. This matrix element is found to be much smaller than the corresponding 
matrix element in the color-singlet model, which measures the probability for the formation 
of i/j from a pointlike cc pair in a color-singlet 3 Si state. In spite of the fact that the 
color-octet matrix element is intrinsically smaller, the color-octet term can dominate over 
the color-singlet term in the gluon fragmentation function because the color-singlet term is 
suppressed by a short-distance factor of a 2 s {m c ). 

In order to establish that the color-octet mechanism plays a significant role in the produc- 
tion of charmonium, it is necessary to find other production processes in which it dominates 
over the conventional color-singlet mechanisms. In this Letter, we point out that there is 
a dramatic signal for color-octet processes in inclusive ip production in e + e~ annihilation. 
The color-octet terms are suppressed relative to the color-singlet term over most regions of 
phase space, but they are enhanced near the upper endpoint of the ip energy distribution. 
The distribution of the angle 9 between the 3-momentum of the ip and the beam axis has 
the form 



where E is the energy of the if) in the center of momentum frame. The shape of the angular 
distribution is characterized by the parameter A(E), which varies with E. The color-singlet 
model predicts that A(E) should approach a value close to —1 near the upper endpoint of 
E. We show in this paper that, if color-octet contributions dominate near the endpoint, 
A(E) will change sign and approach a value close to +1 at the endpoint. This change in the 
sign of A would provide a dramatic signature for color-octet production mechanisms. 



form ip |J. This term corresponds to the gluon forming a cc pair in a color-octet 3 Si state at 




(1) 
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Using the formalism of Ref . , the inclusive cross section for the production of a char- 
monium state H with momentum P can be written in a factored form: 



da(e + e~ -> H(P) + X) = ^rfa(e + e 



cc(P,n)+X) (<%), 



(2) 



rt 



where da is the inclusive cross section for producing a cc pair with total momentum P, 
vanishing relative momentum, and in a color and angular-momentum state labelled by n. 
This parton cross section is a short- distance quantity that can be calculated as a perturbation 
series in a s (m c ). The long-distance factor (O^) is the matrix element of a local 4-fermion 
operator in nonrelativistic QCD which represents the inclusive probability for the formation 
of the meson H from the pointlike cc state n. All dependence on the particular quarkonium 
state H resides in the matrix elements. The relative magnitudes of the matrix elements 
(0%) can be estimated using scaling rules that predict how they scale with the quark mass 
m c and with the relative velocity v of the quark. 

In the case of P-wave states, there is a color-octet matrix element that contributes at the 
same order in v 2 as the leading color-singlet matrix element. The color-octet term is needed 
for perturbative consistency, because it cancels an infrared divergence in the color-singlet 
term |§j . In the case of S-waves, the leading color-octet matrix element is suppressed by v 4 
relative to the leading color-singlet matrix element. It is therefore perturbatively consistent 
to neglect the color-octet term. However, the suppression factor v 4 is only about 1/10 for 
charmonium, so the color-octet term is not necessarily negligible numerically. 

The predictions of the color-singlet model for the inclusive production of if) in e + e~ 
annihilation have been explored thoroughly P~ PH - in this model, the sum (0) is truncated 
to a single term whose matrix element is proportional to |P(0)| 2 , where R(0) is the radial 



cc pair is produced in a color-singlet 3 Si state with vanishing relative momentum. The 



over the energies of the recoiling gluons. The complete expression for the differential cross 
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wavefunction at the origin. In the notation of Ref. , the matrix element 
(9/27r)|#(0) | 2 . At leading order in a s , the parton process is e + e — ► c 




ccgg, where the 



corresponding differential cross section da[ 1 can be obtained from Ref. [ITJ by integrating 



section is too lengthy to reproduce here, but it becomes relatively simple at the endpoints 
of the energy spectrum. Near the lower endpoint E min = 2m c , the distribution is isotropic, 
corresponding to A = 0. The upper endpoint is -E max = (s + 4m^)/(2y / s), where s is the 
square of the center of mass energy. At this endpoint, the distribution approaches 



where r = 4m 2 /s. For comparison with data, we can identify 2m c with the mass of 
the if). Thus the color-singlet model predicts that A(E) should vary from at E = M$ to 
-{s - Ml)/{s + Ml) at E mSLX = (s + Ml)/(2s/s). The prediction for A(E) as a function of 



the if) momentum = J E 2 — M 2 at \/s = 10.6 GeV is shown as a solid line in Fig. 1. The 

predicted value at the endpoint is A = —0.84. 

Color-octet contributions to the differential cross section (0) are suppressed by a factor of 
v 4 from the matrix element (0%). They can be competitive with the color-singlet terms only 
if there is some enhancement from the short- distance cross section da. Such an enhancement 
does occur in the endpoint region. At leading order in a s , the color-octet contribution comes 
from the parton process e + e~ — > ccg, which produces cc pairs only in the endpoint region 
and with angular momentum 1 S , 3 Pj, 1 D , etc. From the velocity-counting rules of Ref. [Q, 
the dominant matrix elements at this order in a s are (Of ( 1 5'o)) and (Of ( 3 Pj)), J = 0, 1, 2, 
which are suppressed by v A relative to the color-singlet matrix element (Of ( 3 5i)) . The 
matrix element (Of ( 1 5'o)) is proportional to the probability for the formation of a ip from a 
pointlike cc pair that is produced in a color-octet 1 So state, and similarly for (Of ( 3 -Pj)). The 
suppression by v 4 arises because the formation of the color-singlet bound state can proceed 
through a magnetic-dipole transition in the 1 S case and through a double electric-dipole 
transition in the 3 Pj case. The three P-wave matrix elements are not independent at leading 
order in v 2 , but are related by heavy-quark spin symmetry: 



The short-distance cross sections corresponding to these matrix elements can be calculated 






(Ot( 3 Pj)) « (2J+l)(^( 3 P )). 



(4) 



from the Feynman diagrams for e + e 



ccg: 
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dai s °' > „. , IGtc 2 a 2 a s , 

dE^Te = 6{E ~ " + cos ^ ' (5) 

da { s Po) MF-F ^ 327^ 2 a 2 a g (l-3r) 2 (l + cos 2 g) 
ciEdcosfl 1 maxJ 27m3 s 2 2 (l-r) ' 1 J 

^f Pl) A , P _p 327r 2 a 2 a,(l + 2r) + (l-2r)cos 2 g 
dEdcos9~ [ maxj 27m3 s 2 (i _ r ) " ^ 

rf?f P2) 32vr 2 a 2 a s (1 + 6r + 6r 2 ) + (1 - 6r + 6r 2 ) cos 2 

dEdcos9~ { max) 21m\s 2 5(1^7) ' ( ) 

The delta functions should be interpreted as narrow distributions in the energy. In the 
quarkonium rest frame, the width of the delta distribution is of order m c v 2 . This is the order 
of magnitude of the splittings between the radial and orbital-angular-momentum energy 
levels of charmonium, which are roughly 500 MeV. The effect of boosting to the e + e~ center 
of momentum frame is to decrease the width of the delta function by the factor 2m c /^/s, 
which at CLEO energies is approximately 0.3. Thus the delta functions in should 
be interpreted as distributions in E whose widths are on the order of 150 MeV. 

We now compare the magnitudes of the color-singlet and color-octet contributions to the 
cross section. Away from the endpoint region, the color-singlet and color-octet parton cross 
sections da are both of order a 2 . The color-octet terms are therefore suppressed by a factor 
of f 4 from the matrix elements. In the endpoint region, the parton cross sections (§|)-(§D 
for the color-octet terms are of order a s and they are enhanced by a factor of 1/v 2 from the 
energy delta functions. The relative magnitudes of the color-octet terms compared to the 
leading color-singlet term is therefore v 2 /a s in the endpoint region. Thus the color-octet 
terms can be competitive in this region. They will dominate over the color-singlet term if 
the 1/otg enhancement overcomes the v 2 suppression. 

If we assume that the color-octet terms dominate in the endpoint region, we obtain a 
simple and dramatic prediction for the if) angular distribution. In the endpoint region, the 
distribution should approach a form that corresponds to some linear combination of the 
differential cross sections dd% So ^ and J2j(2J + l)^5g Pj ^ ■ If the 1 S term dominates, the 
angular distribution parameter A(E) should approach +1. If the 3 Pj terms dominate, then 
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A(E) should approach the value (3 — 10r + 7r 2 )/(3+6r + 7r 2 ). At CLEO energies, this value is 
approximately +0.62. Thus, regardless of the relative sizes of the matrix elements (Of O-So)) 
and (Of ( 3 -Po)), the value of A is predicted to be positive and close to +1. The predictions 
from color-octet dominance of the endpoint region fall in the hatched region in Fig. 1, with 
the upper and lower lines corresponding to dominance by the matrix elements (Of ( 1 5'o)) 
and (Of ( 3 Po))> respectively. This prediction stands in sharp contrast to the prediction of 
the color-singlet model that A(E) should approach —0.85 in the endpoint region. 

Higher order radiative corrections and relativistic corrections will produce a smooth 
transition between the moderate z region where the color-singlet contributions dominate and 
the endpoint region where the color-octet contributions dominate. Relativistic corrections 
to the color-singlet contribution can be taken into account by including additional matrix 
elements in the general factorization formula (0). At order a 2 , the first relativistic correction 
will change the cross section by an amount of order v 2 . This could easily change the color- 
singlet prediction for the asymmetry parameter A at the endpoint by about 30%, but it is 
unlikely to change the sign. At order a^, the relativistic correction diverges logarithmically 
as E — > i? max . The singularity at the endpoint can be factored into one of the color-octet 
matrix elements and it has therefore already been taken into account in our analysis. 

For a more quantitative analysis of A(E), we require phenomenological determinations of 
(Of ( 1 5'o)) and (Of ( 3 Pq)). These matrix elements are independent of the color-octet matrix 
element (Of ( 3 5i)) that plays an important role in gluon fragmentation into if). One process 
which provides information on (Of ( 1 5'o)) and (Of ( 3 -Po)) is photopro duct ion. The cross 
section for ^yN — > if) + X was first calculated in the color-singlet model by Berger and Jones 



12| . The relevant parton process at leading order in a s is 75 — > ccg. As emphasized by 
Berger and Jones, the color-singlet model breaks down in the endpoint region z — > 1, where 
z = E^ J /E 1 is the ratio of the energies of the if) and the photon in the nucleon rest frame. 
From the point of view of the factorization approach of Ref. 0], this breakdown of the color- 
singlet model can be attributed to contributions from color-octet processes in the endpoint 
region z — > 1. At leading order in a s , the relevant parton process is 7^ — > cc, which produces 



color-octet cc pairs in the endpoint region with angular momentum quantum numbers x Sq 
and 3 Pj. Thus it might be possible to use data on photoproduction in the nearly elastic 
region z — > 1 to constrain the matrix elements (Of ( 1 5'o)) and (Of ( 3 P/)}. Unfortunately, the 
endpoint z — 1 for photoproduction is dominated by elastic and diffractive ^'s, which cannot 
be described adequately within our factorization framework. The color-octet mechanism 
produces inelastic ip's with z near 1, and it is not clear whether these can be separated from 
the diffractive contribution. In order to eliminate the diffractive contribution, it may be 
necessary to go to higher momentum transfer and determine the color-octet matrix elements 
from electroproduction data. 

We have shown that the inclusive production cross section for J /if) in e + e~ annihilation 
may be dominated by color-octet production mechanisms near the upper endpoint of the if> 
energy spectrum. The signal for the color-octet contributions is a dramatic change in the 
angular distribution of the if), which has the general form 1 + A(E) cos 2 9. The prediction 
of the color-singlet model is A « —0.84 at the upper endpoint E mSuX . If the endpoint region 
is dominated by color-octet processes, then A(E) should change sign near the endpoint, 
approaching a value A > +0.62. Thus the angular distribution provides a simple and dra- 
matic signature for color-octet production mechanisms. The enormous statistics available 
at CLEO makes it possible to measure the angular distribution as a function of the momen- 
tum of the if). Preliminary results from CLEO indicate that the angular distribution is 
consistent with the color-singlet model at moderate values of p^, but that the asymmetry 
parameter A changes sign at the largest values of p^ |L3| in agreement with the prediction 
from color-octet dominance. Such a result would provide strong evidence that color-octet 
processes play a significant role in charmonium production. 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grant DE-FG02-91-ER40684. 
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FIGURES 



Angular distribution parameter A as a function of the momentum of the for e + e~ 
annihilation at y/s = 10.6 GeV. The solid line is the prediction of the color-singlet model. 
The hatched region is the band of predictions in the endpoint region, assuming that this 
region is dominated by color-octet production mechanisms. 
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